arXiv:1504.04154v2 [gr-qc] 20 Apr 2015 


Interaction between Maxwell field and charged scalar field in de 

Sitter universe 

Cosmin Crucean and Mihaela-Andreea Baloi 
West University of Timi§oara, 

V. Parvan Ave. 4 RO-300223 Timi§oara, Romania 

Abstract 

We study the theory of interaction between charged scalar field and Maxwell field in de Sitter 
background. Solving the equation of interacting fields we define the in-out fields as asymptotic 
free fields and construct the reduction formalism for scalar field. Then we derive the perturbation 
expansion of the scattering operator. The first order transition amplitudes corresponding to particle 
production from de Sitter vacuum and pair production in an external field are analysed. We show 
that all these effects are important only in strong gravitational fields and vanish in the flat limit. 

PACS numbers: 04.62.-|-v 


1 


I. INTRODUCTION 


In de Sitter geometry the problem of particle production was considered by many au¬ 
thors. These studies show that the space expansion generate particle production using 


various methods 


15 


17- 


32l |. However a detailed study of interaction between scalar held 


and Maxwell held minimally coupled with gravitation using perturbative methods and exact 
calculations received little attention. We focus our attention on this subject because the last 
astronomical observations shows that the expansion of the Universe is accelerating, thus in¬ 
creasing the interest on de Sitter model. The study of interactions between bosons and scalar 
helds on a curved background will help us to understand how the scalar electrodynamics 
from hat space will be modihed when a gravitational held is present. Important results such 
as normalised solutions was obtained for Dirac held js], 12| and scalar held jl^, 23|,[w, 11 1 
on de Sitter geometry. In the case of the electromagnetic held the situation is simple if one 
recalls the conformal invariance between de Sitter and Minkowski metrics. This means that 
the theory of free electroma gnet ic held is the same as in Minkowski geometry if we use the 


metric with conformal time 


14| . In the present paper our we will study on the interaction 


between charged scalar held and electromagnetic held. 

We want to study the ehect of electromagnetic interaction upon particle production in 
de Sitter expanding universe. The intention of this paper is to develop the theory of the 
scalar QED on de Sitter spacetime and to use this formalism for studying the following 
processes: vacuum (p + (p* in an external Coulomb held and vacuum (p + (p* + ■y. Here 
ip represents the charged scalar particle, ip* represents the charged scalar antiparticle and 7 
the photon. The annihilation of scalar pair into vacuum in the presence of Coulomb held 
and annihilation of the triplet scalar pair and a photon into de Sitter vacuum represents the 
time reversed processes. Our results prove that these probabilities are nonvanishing only in 
strong gravitational helds. Our exact calculations allows us to study the behaviour of our 
amplitudes/probabilities for large/small values of the expansion factor, making in this way 
the connection with the scalar QED from Minkowski space. We must specify that in our 
calculations we will use the solutions of Klein-Gordon and Maxwell equations which have 
a dehned momentum. These reactions are forbidden in hat space theory because of the 
energy-momentum conservation laws. In de Sitter case the energy is no longer conserved 
and we must expect to obtain a nonvanishing probability for these processes. 
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We start in the second section with a review of the results concerning the free scalar field 
theory on de Sitter space. In the third section we construct the theory of interaction between 
charged scalar held and electromagnetic held on de Sitter space and then our attention will be 
focused on developing the reduction formalism for scalar held. In this section is also studied 
the perturbation theory which help us to obtain the dehnition of scattering amplitudes. The 
hrst order processes that generate particle production will be studied in sections four and 
hve. Further we study the probabilities that generate particle production using a graphical 
method. Our conclusions are summarized in section six, and in the Appendix we give the 
main steps that help us to calculate the transition amplitudes. 


II. FREE SCALAR FIELD 

The exact solutions of the free Klein-Gordon equation on de Sitter space-time^ was ob¬ 
tained in 23|. The calculations was done using the de Sitter line element 

1 


flfl 


ds^ = df - = 


{Ojtcf 


(^dtl — dx‘^) , 


( 1 ) 


where u is the expansion 


Gordon equation reads jl3| : 


actor. In the chart {t, x\ with Gartesian coordinates the Klein- 


— e -I- ?>udt + m?) (j){x) = 0 . 


( 2 ) 


It is known that the Klein-Gordon equation ([2]) on de Sitter space can be analytically solved 


in terms of Bessel functions 


10U23ll . In our study we will use the normalized solutions of 


positive frequencies obtained in jl3| : 

where H^{z) is the Hankel function of hrst kind, p = |p| and the notations are: 

/ m 

( 4 ) 


k=Ju^-j, 


m 

u = — 


with m > 3u/2. The fundamental solutions of negative frequencies are f^{x). For t —oo, 
the modes dehned above behave as a positive/negative frequency modes with respect to the 
conformal time G = 


ff{x) ~ e fi{x) ~ e 


iptc 


( 5 ) 
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This is the behavior in the infinite past for positive/negative frequency modes which defines 
the Bunch-Davies vacuum. 


The fundamental solutions of Klein-Gordon equation, satisfy the ortonormalizations re¬ 
lations Q: 


i j d^x{-gY/^ f^{x) dt = 

-i j d^x{-gy/‘^ ff{x) dt = ^Hp-p'), ( 6 ) 

i J d^xi-gY^^ ff;{x) dt fpYx) = 0 , (7) 

where the integration extends on an arbitrary hypersurface t = const and {—gY^"^ = 

The completeness condition satisfied by these solutions reads [l3|: 

i j d^p fp{t,x) dt fp{t,x') = e~^‘^^S^{x - x'). (8) 


The quantization can be done considering the plane wave in momentum representation: 

(/)(x) = 0(+)(a:) + (/>(-Yx) = j d^p [fp{x)a{p) + fp{x)b^{p)] , (9) 

where the particle (a, a^) and antiparticle {b, b^ operators satisfy the standard commutation 
relations in momentum representation: 

[a{p), aYp')] = [b{p), bYp')] = S^p- p') ■ (10) 

These operators can be calculated using the inversion formulas: 

a{p) = i J d^x{-gY^‘^fpix) dt 4>{x) 

b{p) = i j d^x{-gY^‘^fp{x) dt (f)Yx). (11) 

All these operators act on the Fock space supposed to have an unique vacuum state |0) 
accomplishing 

a(p)|0) = 6(p)|0) = 0 , (0|a^(p) = (0|6^(p) = 0 . (12) 

In this way the field 0 is correctly quantized according to the canonical rule |^: 

[(l){t,x),dt(j)Yt,x')] = ie~^^6^{x — x'). (13) 
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In the quantum theory of helds it is important to study the Green functions related to the 
partial commutator functions (of positive or negative frequencies) dehned as: 

D^^\x^x') = (14) 

and the total one, D = + D'^~\ These functions are solutions of the Klein-Gordon 

equation in both sets of variables and obey [D^^\x, x')]* = D^^\x,x') . Then the partial 
commutator functions will have the expressions: 

D^+\x,x')=i j £p fp{x)fp{x'), 

D^-\x,x') = -i J £pfp{x)fp{x'). (15) 

These functions help us to introduce the principal Green functions. In general, G{x, x') = 
G{t,t',x — x') is a Green function of the Klein-Gordon equation if satishes |^: 

{pi — + 3udt + mP G{x, x') = e~^‘^^6^{x — x'). (16) 

The properties of the commutator functions allow us to construct the Green function just 
as in the scalar theory on Minkowski spacetime. We assume that the retarded, D^, and 
advanced. Da-, Green functions read: 

Dji{t,t',x — x') = 9{t — t')D{t,t', X — X'), (17) 

DA{t,t',x — x') = — 9{t'— t)D{t,t', X — X'), (18) 

while the Feynman propagator, 

DF{t,t',x — x') = i(O|T[0(x)0^(a;')] |0) 

= 9{t — ,x — x') — 9{t'— t)D^~\t,t',x — X'), (19) 

is dehned as a causal Green function. Thus we have now all the elements that we need for 
begin our study related to the reduction formulas for scalar held on de Sitter space-time. 

III. SCALAR QED 

The goal of the quantum held theory is to describe the dynamics of interacting particles. 
We are interested to obtain the amplitudes for single particles to propagate in a curved 
space-time and the transition amplitudes for interacting particles between diherent initial 
and hnal states. Our main attention in this paper will be focused on interaction between 
charged massive scalar held and electromagnetic held on de Sitter space-time. 
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A. Interacting fields 


In this work we use for our calculation the same formalism as in jsl-js], and we adopt 
the minimal coupling corresponding to the classical interaction of a point charge as the 
prescription for introducing electro dynamic couplings. Also we will note the interacting 
fields with 0(a;) and A^{x). These fields will satisfy a system of equations that can be 
obtained from the following Lagrangian density: 

A] = j d^x[Csc{ 4 >) + C,m{A) + CM. A)] 

-ie[{dM)(f> - , ( 20 ) 

where the density of lagrangian of interaction reads Cj = 

yJ—ge^(l)^(f)AM^- This relation is obtained after we make the substitution 9^0 —)■ d^(j)—ieA^(l) 
in the lagrangian of free scalar field. From action (1201) . we obtain the system of equations 
for interacting fields on a curved background as follows: 

^ d, (= -^e[(a>t)0 _ + 2eVV2l", 


Vi-9) 
1 


V^ 


df, [VM9^''dy(j) - + e^cpAM. (21) 


The problem of finding exact solutions for the system (|2T]) . which contain coupled nonlinear 
equations for interacting fields proves to be as in the fiat space case, a complicated task. It 
is immediate from here that for finding solutions of the conned equations, we can try to 
follow the same approximation method as in fiat space case ^-js]. The system of coupled 
equations is replaced with one system of integral equations that contain information about 
initial conditions. For that we select the Green function G{x,y) corresponding to one initial 
condition, which help us to write the solution of the second equation from system fl2T|) as 
follows: 


(j){x) = M) + e / d‘^yV^G{x,y) 


V[V^Vy)A^{y)] 


V^ 

+i[Wy)]A^iy) + e(j){y)AM)A^{y)"^, ( 22 ) 

where 4>{x) are free fields. One can verify that (|2^ is the exact solution of the second 
equation from fl^ . applying {Ekg+^‘^) to (f>{x), with the observation that [Ekg+^'^W) — 
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0. We introduce the following convenable notation Ekg = :^ 7 = g^'^dy). Using the 

fact that Green functions must satisfy equation ffT6|) . it is not hard to verify that the solution 
is justihed. 

Since our main scope in this paper is to study transition amplitudes, we want to construct 
the states which give a simple description of the physical system at the time t —oo. 
At this time the particles have not yet interacted with each other and propagate under 
the influence of their self-interactions only. This implies that for constructing a theory of 
interactions, we must seek for operators that create independent particle states with each 
particle propagating with its physical mass. Let us start to construct these operators, using 
which offers us the possibility of constructing free helds, which are asymptotic equal (at 
t —)■ ±oo) with solutions of system lET]) . As we know the retarded Green functions Dji{x, y) 
vanishes at t —)■ —cxo while the advanced one DA{x,y) vanishes for t ^ oo. Thus equation 
fl2^ can be expressed with retarded and advanced functions as follows: 

(/)(x) = ^R/A{x) + e j d%y/^Dii/A{x,y) 

+i[di,(j){y)]A^{y) + e(j){y)Af,{y)A^^{y)y (23) 

Now we can dehne the free helds, (j)R/A{x) that satisfy: 




lim - ^r/a{x)) = 0. (24) 

The free helds (pR and (pA have mass m and are equal at t ± cxo with exact solutions of the 
coupled equations and represent the helds before and after the interaction. It is known that 
the mass of the interacting helds may diher from the mass of free helds. For obtaining the 
mass correction we will ad to the lagrangian of interaction a term of the form —6m‘^(p^<p, 
which help us to rewrite the second equation from system fl?I|) : 




g^'^diyCp — iey/^cpA^] -f m^cp = ie{dfj,(p)A^ + e^cpA^A^ — 6rri^(p, (25) 


where the current j = ie{d^(p)A^^ + e^(pA^A^ — 5m?(p + --^d^{y/—g(pA^) is now the source 
giving rise to the scalar helds. 

n n 

Further we can follow the method from Minkowski case |3|-[8|, where the free helds (pR 
and <pA are dehned up to a normalization constant noted with ^/z. This allows us to dehne 
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the in/ out fields in order to construct the states of single particles, as follows: 


y/zcpinioutix) = (j){x)-e / d^yy/^DR/A{x,y) 






+i[d^(p{y)]A^{y) + e(j){y)Af,{y)A^{y) - 6m‘^(j){y)j. (26) 

The in/ out free fields defined above satisfy Klein-Gordon equation, are written with the help 
of creation and annihilation operators and in addition satisfy the conditions: 


lim (0(x) - y/z(j)in/out{x)) = 0. 

It is useful to express these fields with the help of equation fl25l) . in the form: 
•/^^>m/aui{x) = ^{x) - / d‘‘yX^DR/A(x, !/)|£^*- g (!/) + 


(27) 


(28) 


From the above normalization we can define the creation and annihilations operators as in 
Section II , 

)m/ont ^ ^ ^ 4^in/out{x') 

i,P') in / out ^ ^ ^ 4^ in / outix') • (^9) 

The creation and annihilation operators defined above satisfy the commutation relations 
flT0|) and from that it follows that all the properties of free fields will be preserved. 

Before starting our calculations we make a few remarks about the scattering operator. 
Denoting the vacuum state |0) (which is suppose to be unique), the one particle/antiparticle 
states can be written: 

a\p)in/outYAi = \in/out,l{p)) 

b\p)in/out\^) = \in/out,l{p)). (30) 

If we consider two states |m, a) and |oMt,/3), then the probability of transition from state 
a to state (3 is defined as the scalar product of the two states: {out, (5\in,Q). These are 
just the elements of matrix for scattering operator Spa, which are of physical interest. This 
operator assures the stability of the vacuum state and one particle state, and in addition 
transform any out field in the equivalent in field. 

Our remaining task is to construct and study, general matrix elements which describe the 
dynamical behavior of interacting particles. As in Minkowski theory S-0 , we can construct 
n—particle in and out states as in the free scalar field theory by repeated application to the 
vacuum of a^{p)in/out and b^{p)in/out- 





B. Reduction formalism 


The reduction formulas for scalar particles can be obtained as in Minkowski theory 
8|. Let us consider the amplitude for the process in which particles from in states denoted 
with a, together with a scalar particle l(p) pass in out state, in which we have a scalar 
particle l{p') and particles denoted by (3. The reduction from out state can be done using 
the following relations: 


O'outip ) ) i I d X€ dt (jpoutipd) 0m(^))) 


(31) 


4>mt{x) - 4>in(x) ^ ^ j d'ysf^D(x,y)[EKa{y) (32) 

where we noted D{x,y) = Dr^x^h) — DA_{x,y). The hnal result for reduction of scalar 
particle from out state is: 


{out(3,l{'p')\ina,l{p)) = 6^{p — p'){out(3\ina) 

x{putl3\(l){y)\ina,l{p))d3‘y. (33) 


This method can be also applied for the reduction of particles from in state and the hnal 
result read: 


{out(5,l{p')\ina,l{p)) = 5^{p — p'){out(3\ina) 

x[EKG{.y)+m^]fp{.y)dS. (34) 

Let us resume the above results related to reduction method. We showed that all particles, 
will be replaced by formulas that were obtained in the reduction of one particle. When 
more particles are reduced, in the matrix element appears the time ordered products of 
corresponding held operators. Every particle will be replaced after we apply the above 
reduction method with expressions which depend on the held operator 0(a:). As in the hat 
space case, we can associate a Green function to any process of interaction that must be 
calculated using perturbation theory. 
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C. Perturbation theory 


Expanding the transition amplitndes and matrix elements in a power series in the interac¬ 
tion strength, we can obtain the rnles for compnting the terms in the expansion. The Green 
fnnctions of the interacting helds cannot be calcnlated exact and for that reason we must 
use perturbation methods. The form of the amplitudes obtained from reduction formalism 
allows us to use perturbation calculus as in Minkowski theory S-B- 

For calculating the Green functions it will be helpful to write them in terms of free helds 
as follows: 

(O|T[0(?/i)0(l/2)...i(i/„), S]|0) 


G'(l/i,2/2,-.,2/n) = (O|T[0(?/i)0(?/2)...kl(|/„)]|O) = (0|S|0) 

Then the entire perturbation calculus is based on expansion of the operator S: 


(36) 


oo / ’\n p 

S = Te ~^^= 1 + / 

n=l 


T [Cl {Xi)...Cl {Xn)]d‘^Xi. ..d^Xn, 


(36) 


where the density lagrangian of interaction is (we neglect the second term in the density 
of lagrangian of interaction); Ci{x) = —ie^—g : [{d^(fC)(f) — (9^0):. Each term from 
(1361) corresponds to an order from perturbation theory. Replacing fl36ll in the expression of 
Green function (|35ll we obtain perturbation series which allows us to calculate the transition 
amplitude in any order. 

With this method, the dehnition of transition amplitudes in any order of perturbation 
theory can be obtained. The Green functions from equation fl35l) can be expressed in terms of 
Feynmann propagators considering all the T contractions from Wick theorem. For example 
in the hrst order of perturbation theory the transition amplitude reads: 


((p')|Si|(p)) = -e 




fp'ix) fpi^) 


A^{x)d‘^x, 


(37) 


where the bilateral derivative acts as: g d h = {dh)g — {dg)h. The above expression 
is just the scattering amplitude in the hrst order of the perturbation theory on de Sitter 
expanding universe. With the above method were established the general rules of calculation 
for deducing the expression for the scattering amplitudes in any order of the perturbation 
theory. The results obtained in the last two sections will be applied to calculate the processes 
which generate particle production. 
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IV. PARTICLE PRODUCTION IN COULOMB FIELD 


In this section we want to stndy the process vacuum —)■ ip+ip*, in the presence of Coulomb 
held, using the scalar QED formalism developed in the previous sections. More precisely, we 
want to show that there are nonvanishing probabilities for scalar pair production in the case 
of expanding de Sitter universe. It is well known that the Coulomb held could not generate 
particle production in scalar QED from Minkowski space as a perturbative phenomenon. 
Contrary to this, in the de Sitter case we must expect that this process to be allowed by the 


fact that energy is no longer conserved. The studies done before with fermions indeed s’ 
that there are nonvanishing probabilities for pair generation in external Coulomb held 


now 


16|. 


The ingredients of our calculations will be the solutions of Klein-Gordon equation in the 
momentum basis and the form of Coulomb held in de Sitter geometry. Using the conformal 
invariance we can establish the expression for Coulomb held in this geometry in the natural 
frame Q,Q: 

Zp 

A* = 0. (38) 


A\x) = 


The amplitude corresponding to scalar pair production in Coulomb held on de Sitter 
spacetime, reads: 


{(p)-.(p"')+|Si|0) =-e / \/-g(x) \fsf(x) d, f}(x)] A!‘(x)d*x, 


(39) 


where the notations (p)- stands for a scalar particle of momentum p and negative charge 
and (p')+ denotes the antiparticle with positive charge and momentum p'. 


A. Probability of pair production 

In order to obtain the expression of the probability of pair production, we hrst need to 
calculate the amplitude of transition fl3^ . We use in our calculations the solutions of Klein- 
Gordon equation in momentum basis ([3]). These solutions are written in terms of Hankel 
functions with imaginary index, provided that m/u: >3/2. The spatial integral is the same 
as in the hat case and the result is incorporated in the amplitude formula. The temporal 
integrals can be expressed with the solutions of the free scalar held ([3]). After we pass to the 
new variable of integration s = /u and calculate the bilateral derivative, the temporal 
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integrals become: 


dss^ 


,u{2) 


-p'H\lk{.Ps)- 


d{p's) 


+ pH^^lip's) 


d{ps) 


(40) 


Using the relation for derivative of Hankel functions flFr|) . and the relation between Hankel 
functions and Bessel K functions fl58l) from Appendix, we hnally obtain the expression of 
transition amplitude in the form: 

((p)_,(p'')+is.io) = -^d^ 

dTT'^lp + p'p 

X j-p' y dss^ [K_ik{ips)K_ik-i{ip's) + K_ik{ips)K_ik+i{ip's)] 

+p j dss‘^ [K_ik{ip's)K_ik-i{ips) + K_ik{ip's)K_ik+i{ips)]'^. (41) 

The results of these integrals are discussed in Appendix dSHD- The hnal result for the 
amplitude is expressed in terms of hypergeometric Gauss functions, gamma Euler functions 
and unit step functions as follows: 


((p)-,(p')+|Si|0) = 


Ze^m 


Stt^Ip + p'P \_p 


9{p-p') ( -hk { — ] +gk 


P 


p, 


+ +hk[-^, 


p' 


P 


p, 


(42) 


where the newly introduced functions Pk , hk {^) have the following expressions (with 


the specihcation that pk , hk ^y'j is obtained when p ^ p'): 


9k 


-1 


m\ 
u ) 


/m\ 


P 


pi 


—ik 


r(l — ik)T{2 + ik) 2 F 1 ( 2,1 — i/c; 3; 1 — 


—ik 


r(2 — -ifc)r(l + ik) 2 F 1 I 1, 2 — zfc; 3; 1 — 


-1 




—ik 


r(3)r(i - ik)r{ik) 2 F 1 \ 2,i-ik-i-ik-[j^ 


+ (^}j r(3)r(2 + ik)T{-ik) 2 Fi ^ 1 , 2 + 1 + ik- 

r(3)r(2 - ik)T{ik) 2 Fi ^1, 2-ik]l- ik] 
+ r(3)r(i + ik)T{-ik) 2 F 1 U,i + ik]i + ik] 


—ik 


(43) 
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hk 


m\ 
u ) 


P 


P 

—ik+2 


—ik 


r(l — ik)T{2 + ik) 2 -F 1 ^1,1 — i/c; 3; 1 — 
r(2 - tk)r{i + tk) 2 F 1 U,2-tk-3-i- 



-1 

u 


—ik 


I r(3)r(l — ifc)r(l + ik) 2 F 1 ( 1,1 — ik] —ik] 


r(3)r(2 + ik)r{—l — ik) 2 -Fi ^2, 2 + ik]2 + ik] 
r(3)r(2 — z/c)r(—1 + ik) 2 F 1 ^ 2 , 2 — ik]2 — ik] 
r(3)r(l - ik)T{l + ik) 2 F 1 f 1,1 + ik] ik] 





ik 


(44) 


The result of the integrals in the case p' = p, give terms proportional with delta Dirac 
function S{p — p'). However when we sum all the contributions of each integral with Hankel 
functions these terms cancel between them. 

From equations flT3l) and (jS]), we can observe that the momentum is not conserved in 
the process of scalar pair production in external Coulomb held. The breaking of momentum 
conservation law is due to the external Coulomb held, because the geometry ([T]) has spatial 
translation invariance as an exact symmetry. If one had just the de Sitter background, 
without the Coulomb held, there would still be certain types of particle production as we 
will show in the next section, but no violation of momentum conservation law 16[ |. 

It is clear that the contribution specihc to the de Sitter geometry in our amplitude is 
contained in the functions gk, hk- These quantities depend of the parameter u = m/cu, and 
contain the ehect of space expansion on the pair production process. These functions will 
dehne the probability of pair production, which is the outcome of our perturbative approach. 
The probability for pair production in Coulomb held is obtained by squaring the amplitude: 


n^/ = l((p)-,(p')+|Si|o)|^ = 


bdvr^lp + p'\ 


XI -^Gi.P-p') 


p-i 


+zii^ip'-p) 


P 


p 

9k 1 — 

P 


h. ( 4 

pf 


+ 


+ 


h,\^- 

P 

p 


pi 


P 


p 


P 


-9k - hi ^ -hk ^ gl ^ 


pi 


P 


pi 


P 


45) 


P 


pi 


P 


p 


9Ab]K[j-,]-hk[^A9l\^ 


P 


pi 


P 


P 


The total probability is obtained by integrating after the hnal momenta p, p' the probability 
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given above. This is a complicated task and we restrict to study the properties of the 
probability given in fH6|) . Since our calculations are exact we can study the results for 
probability in terms of the ratio mass of the particle/expansion factor. 

B. Graphical results 

In what follows we plot the probability as function of the parameter u = m/uj for different 
values of the momenta ratio p'/p. These results are presented in the Figs.(IT]12]). 



FIG. 1: F as a function of u for p/p' = 0.1. 

Our results from the above graphs shows the dependence of the probability in terms of 
parameter u = m/o;, for different values of the ratio p/p'. This result was obtained with the 
exact solution of Klein-Gordon equation expressed with Hankel function of imaginary index 
hI^\pz), where k = \/{m/uY — 9/4, solutions that are valid for m/o; > 3/2. In the above 
graphs the parameter u was taken in the interval u G (3/2, cxo). Only for these values our 
analytical and graphical results are valid. We observe that the probabilities have a variation 
in terms of the ratio p/p'. For the momenta ratio closer to one the probability is very small 
as shown in Fig.([T]). As the momenta ratio become very small the probability increase. This 
result show that the transitions from vacuum in which the momenta ratio are very small will 
be favoured in external Coulomb held (see Fig. ([3])). This is equivalent to say that there are 
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FIG. 2: P as a function of u for p/p' = 0.001. 



FIG. 3: P as a function of u for p/p' = 0.00001. 

bigger probabilities to produce pairs in which the momenta of the particle p, is very small 
comparatively with the momenta of the antiparticle p' {p' » p). 

Another observation that emerges from graphs Figs.([T]l3]) is that the probability decreases 
as the ratio m/u increase. In other words, as the expansion parameter become smaller the 
probabilities for pair production reduces to zero. When the parameter u = oo, which 
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corresponds to the flat space limit (zero expansion factor), the probabilities are zero as we 
can observe from Figs.([I]l3]). This result is also obtained when we replace the ratio m/cj = oo 
in the expressions for amplitudes and probabilities. We can draw the conclusion that in the 
Minkowski limit there are no production of scalar particles in external Coulomb held. 

Further we can address the problem of separation between particles and antiparticles. 


The probability is proportional with the factor 
expanded as: 




-. The modulus from nominator can be 


\p + p'\'^ = {p‘^+p'‘^ + 2pp'cosf3Y, (46) 

where (3 is the angle between the two momenta vectors. From here we can see that for /? = 0 
the momenta of particle and antiparticle are parallel and have the same orientation. In this 
case the probability is minim and the pair will annihilate in vacuum. The opposite case 
is when = vr, which corresponds to the situation when the momenta of the particle and 
antiparticle are parallel and have opposite orientation. This is the case when the probability 
is maxim. Only in this case the particles and antiparticles could separate. 


V. ELECTROMAGNETIC PARTICLE PRODUCTION 

The process: vacuum —)■ <^9 + (p* + 7 , in which scalar pair and a photon is generated from 
de Sitter vacuum is the topic of this section. This process is forbidden in scalar QED from 
flat space by the momentum and energy conservation laws. In de Sitter case the energy is 
no longer conserved and the probability for production of the triplet is no longer zero. 


A. Transition probability 


For computing the transition amplitude we will use the solution of the Klein-Gordon 
equation and the solution of Maxwell equation on de Sitter background, written in the 
momentum basis. The theory of Maxwell held in de Sitter geometry is simple is we recall 
the conformal invariance (where is the vector potential in flat space). 

Then the fundamental solution of Maxwell equation in de Sitter geometry in momentum 
basis is [l^: 


= e 


—2uj t 


(27r)3/2 ^ 




(47) 
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The transition amplitude for this process reads: 


Up)-, ip')+, (<?, A)|Si|0) = -e f ^/-g{x) \fp,{x) di /|(x)] w\x)*d‘^x. (48) 


Our calculations will be valid only in the case m/u >3/2 and we will use the solutions 
of the Klein-Gordon equation expressed in terms of Hankel functions of imaginary index. 
Replacing all the quantities of interest in our amplitude we obtain that the spatial integrals 
give the delta Dirac function 6^{p + p + q), which assures the momentum conservation in 
this process. In the case of the temporal integrals we pass from Hankel functions to Bessel 
J functions, arriving in this way at integrals fl631l64p . which are discussed in Appendix. The 
transition amplitude can be written passing to the new variable of integration s = 


((p'),(p),(?,A)|Si|o) = - 
x[e'^^J_ik{ps)J-ik{p's) + e 


ie7i{p-p')el{q) ^3 

4i/2g (27r)^A sinh^(7rA;) 

~'"^Jik{ps)Jik{p's) - J-ik 


{p + p +q) / dsse 

Jo 

{ps)Jik{p's) - Jik{ps)J-ik(p'p)]- 


(49) 


The final result is expressed in terms of Gauss hypergeometric functions 2 F 1 and Appell 
hypergeometric functions of double argument F^: 


((p'),(p),(<f,A)|Si|o) = 


ZCTl 


-Xp- p')el{q)S^{p + p' + ^ 


4x/2(27r)3/2 

x[hk{p,p',q) + h-k{p,p',q) + fk{p,p',q) + f-k{p,p',q)]- 


(50) 


The functions which dehne our amplitude are dehned as follows: 

ik 


fk{p,p',q) = 


p) ( 1 .1 1 + a. 1 - kt, ^ + , 0 . ^ + iO I (51) 


1 / 3 v'“ v 

~] “r'^T7"TG572^4 l + ifc;^ + i 0 ,^ + i 0 

p J 7r/csmh(7rfc)g®/^ V ^ q q 


3^ , -n 

p2 2 


hk{p,p',q) = 


ie + 1/4) 

4(pp/)3/2 gjj^h^(7r/c) cosh(7rfc) 
^ .3 „ 1 + z 

+ 2-^1 ( 2 2 ~ —2— 




(52) 


In equation fl51ll52p . the newly introduced notation 2 ; is dehned as: 

p2 _|_ p /2 _ g,2 


z = 


2pp' 


(63) 


17 
















We observe that the momentum is conserved in the process vacuum —)■(/) + 99 * + 7 . In 
the previous section we obtained that the process of pair production in external held broke 
the momentum conservation law. From here we can draw the conclusion that the particle 
production processes in the presence of gravity, preserve the momentum conservation law as 
long as there are no external helds (electric helds or magnetic helds). 

The probability is obtained squaring the amplitude and summing after the photon polar¬ 
izations A. Since our amplitude is proportional with a delta Dirac term , by squaring 

this term we obtain \S^{p)\‘^ = VS^{p). In this way we will dehne the probability per unit 
of volume: 


'Pi^f - nYl 


l^l((p'),(p),(g,A)|Si|0)p 






e‘^[{p — p')e\{q)]‘^ 6^{p + p' + q 
64 {2tiY 

\2 , \ r / „ m 2 


\hk{.P,p\q)? + |h_fe(p,p',g)p 


+ \fk{p,p\q)\ + \ f-k{p>,p\q)\ + hk{p,p',q)h*_k{p,p',q) + hk{p,p',q)fk{p,p',q) 
+hk{p,p',q)flkiP:P':q) + h-k{p,p',q)hl{p,p',q) + h-k{p,p',q)fk{p,p',q) 
+h-k{p,p',q)flkiP:P':q) + fk{p,p',q)hlip,p':q) + fk{p,p',q)h*_k{p,p',q) 

+ fkiP:P’,q)f*-kiP:P':q) + f-kiP:P’:q)KiP:P':q) + f-kiP:P’:q)h*_kiP,P’,q) 

+f-k{p,p',q)fk{p,p',q)]- (54) 


The probability depend on the factor: \{p' ~ P)^*x{q)\‘^- Using that the polarization 

vectors are orthogonal on photon momentum q-e\{q) =0 and for any polarization A = ±1, 
they satisfy: 

e\{q)i £\{qTj = Sij - ^ , (55) 

A ^ 

we can express the sum from probability taking into account that the momentum conserva¬ 
tion as given by the delta Dirac term gives p = —{p' + q) '■ 

\{p’- p)el{q)\‘^ = 5^|2p'e*(g)|2 = 4 

A A 

where 6qpi is the angle between the momenta vectors p' and q. This result shows that the 
probability will increase when the angle between the two vectors approach the value 7r/2 
and vanishes for 6qpi = 0,7i. The conclusion is that the scalar particle with momentum p' 
and the photon will be emitted on different directions. The same conclusion emerges if we 
use the angle between p and q. 



p' — 


{p' ■ q. 


= Ap 


12 ■ 2 
sm I 


\p'i 


(56) 
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The Minkowski limit can be obtained when the expansion factor vanishes or equivalently 
when u = m/u = oo. In this limit the probability for generation of the triplet vanishes, 
thus obtaining the well known result from flat space where this process is forbidden by the 
energy conservation law. 

B. Graphical results 

In this subsection the probability given in Eq. fl54p will be analysed in terms of the pa¬ 
rameter u = m/u, using numerical values for the momenta p ,p' ,q ■ 



FIG. 4: P as a function of u for p = 2, p' = 2, q = 1- 

For plotting the probability we take the Appel hypergeometric functions F 4 to be equal 
with one, since these functions are less studied. Using the definition of functions F 4 , one can 
observe that for large u these functions are very convergent, and from here our approximation 
is well justified. Also one can observe from (l5Tll . that the shape of our graphs will be mainly 
determined by the factor (fc sinh( 7 rfc))“^, which is very convergent at large u. We observe 
from our graphs that the probability drops rapidly to zero when the expansion factor u 
decreases. From graphs Figs.dMS]), it is clear that in the Minkowski limit the probability 
vanishes confirming in this way our analytical results for u = 00 . 
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FIG. 5: P as a function of u for p = 1, p' = 2, g = 2. 

VI. CONCLUSIONS 

We studied in this paper the interaction between Maxwell held and charged scalar held 
in de Sitter geometry using the perturbative methods as in hat space case. We established 
the equations of the interacting helds in the hypothesis of minimally coupling. Using then 
the solutions of the equations for interacting helds we construct the reduction formalism 
for the scalar held in de Sitter geometry. The dehnition of transition amplitudes was also 
obtained using the same perturbative methods as in hat space case. We compare two kinds 
of particle production in de Sitter geometry. Our study show that the particle production 
processes in external helds, do not preserve the momentum conservation law. Contrary to 
this we obtain that the mechanism of production of particles from de Sitter vacuum always 
preserve the momentum conservation law. Our analytical and graphical results prove that 
the probability of particle production is signihcant only in strong gravitational helds and 
vanishes in the hat space limit ca = 0. Indeed a direct perturbative calculation in Minkowski 
scalar QED prove that the amplitudes corresponding to our processes are zero. 
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VII. APPENDIX 


The main steps leading to amplitude of pair production in Coulomb field (1^ . will be 
detailed in this section. First, we will use the relation for the derivative of Hankel functions 


35|: 


_ „( 2 ) 




dz * M+1 

and then the connection between Hankel functions and Bessel K functions 


33 


35|: 




(57) 


(58) 


arriving in this way at integrals of the type 


3^: 


i K,(az)iaiz) = r ^^-jr(^-^- 


2 y V 2 

Re(a + b) > 0 , Re(A) < 1 — jRe(ju) | — \Re{v) \ 


(59) 


In our case A = —2 and the second condition for convergence is satished. We also observe 
that in our case a, b are complex and for solving our integrals we add to a a small real part 
a —)■ a + e, with e > 0 and in the end we take the limit e —)■ 0. This assure the convergence 
of our integral and will correctly define the unit step functions and Qk functions. 

The hnal form of the functions , hk which dehne the amplitude of pair production in 


Cou 


33, 


omb field is obtained if we use the following identity between hypergeometric functions 
35|: 

r(c)r(c-a-6) 


2Fi{a,b]C]z) = 


r(c — a)r(c — b) 

+ b - c) 


2 ^ 1 ( 0 , b] a + b — c + l',l — z) (60) 

2 ^ 1(0 — a, c — b; c — a — b + 1; I — z). 


Tia)m 

For the second process of particle production from de Sitter vacuum, we use the well 


known formula for Hankel functions in terms of Bessel J functions 


33 


35j: 


HlP(z) = 


i sin(7r/i) 


^(2)7 A ^ 

^ ^ ?sin(7rp) 


(61) 
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together with; 


(62) 




K^U)= 


and we arrive at the following integrals 35| : 



dxx^^^Ki {cx)Jy{px)Jy{j)'x) = 


which depends on variable u = 

The second type of integrals used are 


n/2 


d 


\/^(pp')T 2 du 


Qu-d'^) 


35]: 


(63) 


{cx)Jy{ax)J_y{bx)dx = 


a \ ^ smvrz/ 


\/^i> 

, 3 a? lP‘\ 

x ^4 ( 1, 1 + z/, 1 - z/; - —, —^ . 


(64) 


These integrals are convergent for Re{c) > 0. For this reason we must replace c —)■ e — iq , e > 
0, when we evaluate these integrals. In the end we take the limit e ^ 0. 

For establishing the hnal results we use the following relations for Legendre functions of 
second kind 35j]: 


Quiz ± iO) = 


TT 


2sin(7rz/] 


[eT‘-"'P„(2) - F„(-j)] , 


(65) 


when — 1 < z < 1 and their relations with hypergeometric functions 2 F 1 351 


d ^ , , 7rz/(z/ + 1) 

—Quiz±t0) = ———- 

dz 4sm TTZ/ 


2-^1 ( 1 — z^, 2 + z/; 2; —-— ) +2 F\ ( 1 — z/, 2 + z/; 2; 


1 + ;^ 


(66 


The Legendre functions of hrst kind can be written in terms of hypergeometric functions 


2F1 


35|: 


Puiz) = 2 F 1 -Z/, 1 + z/; 1; 


1 -z 


(67) 
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